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The method of truncated Mellin moments in a solving QCD evolution 
equations of the nonsinglet structure functions F™ s (x, Q 2 ) and g^ s (x, Q 2 ) 
is presented. All calculations are performed within double logarithmic ln 2 x 
approximation. An equation for truncated moments which incorporates 
ln 2 x effects is formulated and solved for the unintegrated structure function 
f N (x,Q 2 ). The contribution to the Bjorken sum rule coming from the 
O 1 region of very small x is quantified. Further possible improvement of this 

D , approach is also discussed. 
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. 1. Introduction 

Structure functions play a central role in the perturbative QCD. Ex- 
perimental measurements of the spin dependent and unpolarized structure 
functions of the nucleon allow verification of sum rules and determination 
of free parameters of the input parton distributions. From the other side, 
theoretical analysis within perturbative methods investigates the available 
experimentally region of the variables x and Q 2 and the interesting very 
small x region (still unmeasurable) as well. In this low x region QCD 
predicts a strong growth of structure functions with decreasing x - the lon- 
gitudinal momentum fraction of a hadron carried by a parton. Small x 
behaviour of the singlet unpolarized structure functions, driven by gluons, 
is governed by BFKL 1 or CCFM [3] equations which generate the 
steep x~ x (A ~ 0.3) shape of Fi(x,Q 2 ) @]. For the nonsinglet unpolarized 
FJ? S (x,Q 2 ) the driving term at small x is a nonperturbative contribution 
of the A2 Regge pole F^ s ~ x ' 5 ,41 which dominates even over the a s ln 2 x 
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effects. These double logarithmic terms control however the small x be- 
haviour of both nonsinglet and singlet spin dependent structure function 
gi(x, Q 2 ) 0], Knowledge about structure functions at very low x is very 

important. The sum rules, which can be verified by experiments, concern 

1 l 
moments of structure functions / dxx n ~ 1 gi(x,Q 2 ) and J dxx 1l ~ 1 F2(x,Q 2 ) 

o o 
and hence require knowledge of g\ and Fi over the entire region of x £ (0; 1). 
The lowest limit of x in present experiments is about x ~ 10~ 5 so in the- 
oretical analysis one should extrapolate results to x = and x = 1. More 
important is however the extrapolation to x —* 0, where structure functions 
grow strongly than the extrapolation to x = 1, where structure functions 
are equal to 0. The limit x — ► which implies that the invariant energy W 2 
of the inelastic lepton-hadron scattering becomes infinite (W 2 = Q 2 {\ — 1)) 
will never be attained experimentally. So we will really never know "what 
it happens" with structure functions at x — > 0. This situation is however 
not quite hopeless. One can combine the QCD perturbative analysis in the 
very interesting small x region with experimental data without uncertainty 
from the region where x —* 0. It could be achieved through dealing with 
truncated moments of the structure functions, where one takes an integral 
over xo < x < 1 instead over the whole region < x < 1. In usually 
used method of solving evolution equations in QCD one takes Mellin (full) 
transforms of these equations what gives possibility of analytical solutions. 
Then after inverse Mellin transform (performed numerically) one can ob- 
tain suitable solutions of original equations in x space. In this way e.g. in a 
case of DGLAP approximation, the differentio-integral equations for parton 
distributions q(x, Q 2 ) after Mellin transform change into simple differential 
and diagonalized ones in moment space n. The only problem is knowledge 
of input parametrizations for the whole region < x < 1 what is necessary 
in the determination of moments of distribution functions. Using truncated 
moments approach one can avoid uncertainty from the unmeasurable x — ► 
region and also obtain important theoretical results incorporating perturba- 
tive QCD effects at small x, which could be verified experimentally. Trun- 
cated moments of parton distributions in solving DGLAP equations have 
been presented in [5]. Authors have shown that the evolution equations for 
truncated moments though not diagonal can be solved with good precision. 
This is because each n-th truncated moment couples only with n + j -th 
(j > 0) truncated moments. In our paper we adopt the truncated moments 
method to double logarithmic ln 2 x resummation. However a technique we 
use in this approach is quite different because suitable integral equations 
which correspond to ln 2 x terms are of course different from the DGLAP 
ones. As a result we obtain the equations for truncated moments of the un- 
integrated structure function f (x, Q 2 ), where each n-th moment (n 7^ 0) 
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couples only with itself and with 0-th moment. For fixed coupling constant 
a s the result for n-th truncated moment can be found analytically. The 
purpose of this paper is to start the truncated moments method in the case 
of perturbative QCD formalism, describing double logarithmic ln 2 x terms in 
the nonsinglet structure functions F 2 S and g± s . In the next section we re- 
call the approach which resums the double logarithmic terms. The integral 
equation for the nonsinglet unintegrated quark distributions f N (x,Q 2 ) is 
presented. In point 3 the equation for truncated moments of f NS (x,Q 2 ) 
within ln 2 x approach is derived. This equation is solved analytically for 
fixed a s . Agreement for the limit case xq — > (full moments) of our re- 
sults is shown. Results for truncated moments of F^ s and g^ s for simple 
Regge-type input parametrizations and different xq are presented in section 
4. We calculate also the contribution to the Bjorken sum rule coming from 
the region of very small x. Finally in section 5 we summarize our results 
and discuss further possible improvement of our treatment. 



2. Idea of double logarithmic ln 2 x resummation for the 
nonsinglet unpolarized and polarized structure functions of 

the nucleon. 

It has been noticed [1], [5] that the nonsinglet structure functions of 
the nucleon in the small x region are governed by double logarithmic terms 
i.e. powers of a s ln 2 x at each order of the perturbative expansion. This 
contribution to the ln 2 x resummation comes from the ladder diagram with 
quark and gluon exchanges along the chain - cf. Fig.l. In contrast to 
the singlet spin dependent structure function, for the nonsinglet one the 
" bremsstrahlung" nonladder corrections vanish for the unpolarized struc- 
ture function and are negligible for the spin dependent one [5], [7j. In this 
way we do not need to take into account the nonladder diagrams in the case 
of nonsinglet (both polarized and unpolarized) structure functions. The 
Regge theory |B|, which concerns the Regge limit: x — > predicts singular 
behaviour of nonsinglet, unpolarized distributions and nonsingular (flat) 
shape of nonsinglet polarized ones: 

q NS ~ x-°- 5 (2.1) 

A^-^-sr - 5 (2.2) 

where q NS denotes nonsinglet unpolarized and Aq NS nonsinglet polarized 
quark distributions. The nonsinglet part of the unpolarized structure func- 
tion has the form: 



F 2 NS (x, Q 2 ) = F$(x, Q 2 ) - F?(x, Q 2 ) 



(2.3) 
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Fig. 1. A ladder diagram generating double logarithmic ln 2 (l/a;) terms in the non- 
singlet spin structure functions. 



where 

F 2 (x,Q 2 )= Yl e 2 [xq t (x,Q 2 )+xq t (x,Q 2 )} (2.4) 

i=u,d,s,... 

The spin dependent structure function is: 

g? s (x,Q 2 )=g p 1 (x,Q 2 )-gUx,Q 2 ) (2.5) 

and 

9i(x,Q 2 ) = \ Yl e 2 ^ qi (x,Q 2 ) (2.6) 

i=u,d,s,... 

p and n in above formulae denote proton and neutron respectively, ei is a 
charge of the i-flavour quark. Hence finally we get: 

xq NS = F NS {X}Q 2 ) = ^ {Uml _ dval){XiQ 1) (2.7) 

and 

Aq NS = g? s (x,Q 2 ) = ±(Au val - Ad val )(x,Q 2 ) (2.8) 

where u va i, d va i, Au va i, Ad va i are respectively spin and nonspin valence 
quark distributions in the proton. In the double logarithmic approximation 
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the unintegrated nonsinglet quark distribution function f NS (x, k 2 ) satisfies 
the following integral equation 

1 k 2 /z 

f NS (x, k 2 ) = f NS (x) + a a J ^ J ^rf NS (- z , k' 2 ) (2.9) 

X u2 

where 

* = 17 < 2 - 10 > 

and fo /s (x) is a nonperturbative contribution which has a form: 

fo N = a. j - q NS (z) ~ a s x" - 5 (2.11) 



for nonspin distributions or 

l 

f NS (x) = a s f -Aq NS (z) ~ a. In - (2.12) 

J z x 

x 

for spin dependent distributions. The driving terms of these nonpertur- 
bative contributions q NS and Aq NS are shown in (|2.1I) and (|2.2|) . The 

unintegrated distribution f N (x,k 2 ) is related to the quark distributions 
q ns (A^iVS) via . 

^^^^Sf 5 (2 - 13) 

in the unpolarized case and 

f NS (^) = dg ^ k f ) (2.14) 

in the polarized one. Using a method of the Mellin moment functions one 
can obtain from (|2.9|) the following equation: 



f NS (n,k 2 ) = f ^(n)+- 



n 



2 2 / 2 \ ^ 



/c' 2 k' 2 \ k' 2 



k 2 k 2 



(2.15) 

where the full Mellin moment of a function /(x) is defined as: 



l 

n-l 



Af(n) = / dxx n ~ L f{x) (2.16) 



G 



and in our case: 



f NS (n,k?) = I dxx n - l f N »(x,k z ) 



n-l f NSi 



(2.17) 



As it was shown in [1], equation (|2.15|) for fixed coupling a s can be solved 
analytically. Because we would like later to compare full-moments approach 
with the truncated moments one, let us recall the way to obtain the ana- 
lytical solution of (|2.15|) . A short explanation of this is given in Appendix 
A. Thus we get the solution of (|2.15|) in the form: 



2\T 



(2.18) 



where 



and 



n 
2 



1_(^)2 

n 



n = 2^/a~ s 



(2.19) 



(2.20) 



The inhomogeneous term / (n) in (t2~T5l and (t2~THl) according to (|2~THj) . 
(t2~Tll - (f2~T2l) and behaves as: 



j NS i x 

/o (n) 



OL, 



n{n — 0.5) 
for unpolarized structure functions and 



(2.21) 



Jo N ~ — 

n z 



(2.22) 



for the spin dependent ones. The anomalous dimension of the moment of the 
nonsinglet structure function 7 from (|2.19|) has a (square root) branch point 
singularity at n = rig. This gives the following behaviour of the nonsinglet 
structure functions at small x: 



f NS (x,k 2 ) 



2^ _ x -n 



k_ 



2 \ n o/2 



(2.23) 



and hence also 



2^ _ x -n 



■2 \ n o/ 2 



(2.24) 
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This low x shape ~ x~ n °, where no given in (|2.2()|) is equal to 0.39 for 
d s = 0.038 (a s = 0.18), remains nonleading in the case of the nonsinglet 
unpolarized structure function in comparison to the contribution of the non- 
perturbative A2 Regge pole (j2.1j) . In this way QCD perturbative singularity 
at small x generated by the double logarithmic ln 2 x resummation for non- 
singlet unpolarized quark distributions is hidden behind the leading Regge 
contribution: 

q NS (x,Q 2 ) ~ x~ n ° < ar^ 2 (°) a A2 (0) = 0.5 (2.25) 

Quite different situation occurs for the nonsinglet spin dependent functions, 
where the double logarithmic contribution becomes important: 

Aq NS (x, Q 2 ) ~ x~ n ° > a Al (0) < (2.26) 

This takes place because the nonperturbative Regge part for spin depen- 
dent quark distributions involves a very low intercept a J i 1 (0) < 0. Small x 
behaviour ~ x -^\f^l Q f the- nonsinglet structure functions originating from 
double logarithmic ln 2 x resummation is a very interesting feature. Partic- 
ularly for the polarized structure functions, where the double logarithmic 
analysis enable estimation of parton parametrizations at low x. In the next 
section we introduce the truncated moments method and combine it with 
the ln 2 x approach. This technique will give a novel advantage in the QCD 
perturbative analysis: enable to avoid dealing with the unmeasurable x — ► 
region. 

3. Truncated moments method within double logarithmic ln 2 x 
resummation for the nonsinglet structure functions. 

Truncated moments of parton distributions have been lately used in the 
LO and NLO DGLAP analysis 0. Authors avoid in this way an extrapo- 
lation of well known quark distributions behaviour to the unmeasurable so 
unknown x — > region. Apart from they receive evolution equations for 
the truncated moments, in which n-th truncated moment couples only with 
n + j-th (j > 0) moments and that the series of couplings is convergent, 
what ensures a good accuracy. In our double logarithmic analysis we use 
truncated moments of the unintegrated function f NS (x, k 2 ) and for fixed 
coupling a s we can solve the suitable equation analytically. Let us shortly 
explain this treatment. The truncated n-th Mellin moment of the function 
f NS (x,k 2 ) from (t2~T3l - (t2~T31) is defined as: 

1 

f NS (x ,n,k 2 ) = J dxx^f^&k 2 ) (3.1) 

XQ 



The evolution equation (|2.9j) . generating double logarithmic terms ln 2 x in 
the truncated Mellin moment space takes a form: 



k 2 /x 



f NS (x ,n,k' 2 ) = / JV °(>o,ra) +a 



NS, 



dk 



/2 



k' 2 



x I dyy n - l f NS (y,k' 2 ) I dzz n - 1 ® 



1.2 



k' 2 



J-t) 



1 for t > 
for t < 



where Q(t) is Heaviside's function: 

e(t) = 

and we deal with the x-Bjorken region, where 

x > Xo 

After taking into account the relations: 



(3.2) 



(3.3) 



(3.4) 



dzz^Q^-z 



1 

n 



xo/y 



k 2 



x 



; 

(3.5) 



A; 



k 2 



dzz @ \7j2- z ) =lnJ r + @ ( k k ) ln ^ n = ( 3 - 6 ) 



k' 2 

xo/y x 
one can obtain from QM.2JI : 

ATS, 



X 



k 2 /x 



f NS (x ,0,k 2 ) = f 1 " J (x ,0) + a s [ 

k 2 /x 1 



jW2 W2 



A: 2 



+ 



dft f dy l„ „.fNSf„, iJ2 



k 2 /x 



k' 2 J y 



dk 



12 



\nyf^{yX 2 )-\nx I -Z-f^(x ,0,k' 2 )} 



k' 2 



k 2 



k 2 



(3.7) 
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f NS (x ,n,k 2 ) = f NS (x ,n) + / ^-f NS (x ,n, k' 2 ) 



and 

*\ 

k 2 /xo 2 / \ n k 2 /xQ 

+ J f NS (^' 2 )-* n o J ^f NS Mk'% n^O 

fc 2 fc 2 

(3.8 

For ri/Owe get the following solution (for details see Appendix B): 

/™cW) = /o" s (^,») g) 7 1 + (3-9) 

where 

R = R{n,a s ) = 1 ^ (3.10) 

______ — g 

7 is given m (HH and/o (_o , n) is the inhomogeneous term, independent 
on A: 2 : 

l l 

f NS (x ,n) = f dxx^f^ix) = ^ / — (x n - x%) P0 (x) (3.11) 
J n J x 

x x 

The input parton distribution po(x) in the above formula denotes qQ S (x) for 
the unpolarized case or AqQ S (x) for the polarized one, respectively. From 
(|3.9|) one can read that our solution for the truncated moment f NS (xo,n, k 2 ) 
reduces to (|2.18|) when xq = what must be of course fulfilled. 

4. Some results for truncated moments f NS (xo,n,k 2 ), 
F2 NS (xo, n, fc 2 ), cfi NS (xo,n,k 2 ) in the double logarithmic ln 2 x 

approximation. 

The truncated n-th moment of the unintegrated nonsinglet function 
f N (x, fc 2 ) is given by eqs. (|3.9|) - (|3.11jl . In Figs. 2-3 we plot the moments of 
/ (~, k 2 ) for different n as a function of xq at k 2 = lOGeV 2 . One can see 
that the ratio pf(xo,n) defined as 

^_ f NS (x ,n,k 2 ) 
P/(*o.") = /JV5(Qnife2) (4-1) 
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becomes very large (~ 1) at xq rs 10 4 (pf(xo 



10" 



, n 



1) = 0.997). 



This could be an advice that the truncated moments method is useful in 
the region x > xq ~ 10~ 4 because at lower x the results for the full and 
the truncated moments are practically the same (at least for the Regge in- 
put). From the other side double logarithmic ln 2 x terms from our approach 
become important in the small x region x < 10~ 2 [7j, 0. So taking into ac- 
count both above facts one should choose the limit point xo in the truncated 
moments of f NS , gi S or F^ 3 as 



1(T 2 >x > 10~ 4 



(4.2) 



Truncated moments of the nonsinglet quark distributions are related to the 



Fig. 2. Truncated Mellin moments f NS (xo, n, k 2 ) = J dxx' 1 ^ 1 f NS (x, k 2 ) as a func- 
tion of Xo for different n at k 2 = lOGeV 2 in ln 2 x approach. 



moments of / via: 



1 1 k 2 /x 

n-lJVS/_ ;2 n _ f j__n-l „0NS, 



dxx n - l $ b {x,k z ) = J dxx n - L gr b (x) + 

XO XO 



/•- 

I) 



rlk' 2 



for the spin dependent structure function and via: 

1 1 k 2 /x 

n-l^NSt^ i2\ - f J„„n-2 &0NS 



(4.3) 



dxx n -*F? a (x,k?) = I dxx n -' l F^s{x)+ I ,, n/ 7; k , 2 f NS (x ,n,k' 2 ) 



dk 



rl 



xo 



XO 



1-2 



k ,2 {l + 



(4.4) 
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0.002 0.004 0.006 0.008 



Fig. 3. The truncated Mellin moment f NS (xo,n — l,fc 2 = lOGeV 2 ), similarly as 
in Fig. 2 but in a linear scale. 



for the unpolarized case. This together with (THty - lEHTIl gives the following 
formulae: 



0NS / 



+B(x ,n,k 2 )[g l 0NS (x ,n) - x^ s (x ,0)} 



n- 0NS/ 



I 2 (x ,n,k 2 ) = I dxx^F^ixX) = F2"" J (x ,n) 



0NS, 



+B(x , n + 1, k 2 )[F 2 0NS (x , n) - x n Q +1 F 2 mS \x Q , 0)] 



(4.5) 



(4.6) 



where 



B(xq, n, k 2 



7(gr 



ln-i- 



l + (R-l)x% J 1 + e 



dt- 



=7* 



In -ft 

h 1 



(4.7) 



In Table I we collect results for different n-th moments of F^ (x,k 2 ) and 
gi S (x, k 2 ) functions obtained from (|4.5[l - (|4.6|l . The moments are truncated 
at x = HT 1 , 1CT 2 , 1CT 3 and 1CT 4 . For all moments in Table I k 2 = lOGeV 2 , 
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x 


n 


Ii(x ,n, P) 


h{xQ,n, k 2 ) 




1 


0.144631 


0.031977 




2 


0.039952 


0.010018 


10- 1 


3 


0.014137 


0.003945 




4 


0.006101 


0.001843 




1 


0.219244 


0.038540 




2 


0.043831 


0.010413 


1(T 2 


3 


0.014400 


0.003975 




4 


0.006124 


0.001846 




1 


0.229486 


0.038791 




2 


0.043892 


0.010416 


10~ 3 


3 


0.014402 


0.003975 




4 


0.006125 


0.001847 




1 


0.230692 


0.038801 




2 


0.043894 


0.010416 


1(T 4 


3 


0.014402 


0.003975 




4 


0.006125 


0.001847 



Table 1. Truncated moments Ii(x ,n, k 2 ) = J dxx n 1 g^ s (x,k 2 ) and 
1 

l2(xo,n, k 2 ) = J dxx n ~ 1 F£ ,s (x, k 2 ) for different x and n in the ln 2 x approach. 



a s = 0.18 and 



1 

ONS / \ f 

2 {xo,n) = j 



dxx^F^^x) 



Xo 



(4.8) 



1 

^i 0JVS (x ,n) = j dxx n ^g\ NS (x) 



Xo 



(4.9) 



Input parametrizations g± NS (x) and F 2 0,vs (x) are chosen at kfi = lGeV 2 in 
the simple Regge form: 



g° NS (x) = 0.838(1 - xf 



(4.10) 
(4.11) 
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which satisfy the sea flavour symmetric Gottfried [E2 and Bjorken jllj sum 
rules: 

I 0SR = J^Fi !S (x,k 2 )= 1 - (4.12) 
o 

l 

Ibsr = J dxg? s (x, k 2 ) = ^g A = 0.21 (4.13) 
o 

where qa = 1.257 is the axial vector coupling. In order to determine the 
moment integrals we need knowledge of structure functions over the entire 
region of x. The small x behaviour of structure functions is driven by 
the double logarithmic ln 2 x terms. This ln 2 x approximation is however 
inaccurate in describing the Q 2 evolution for large values of x. Therefore 
double logarithmic ln 2 x approach should be completed by LO DGLAP Q 2 
evolution. Dealing with truncated moments of f NS within unified ln 2 x+hO 
resummation one encounters difficulties because suitable evolution equations 
are of course more complicated than in the pure double logarithmic approach 
here presented. Probably the only successful method to solve the evolution 
equations for truncated moments in Zn 2 x+LO treatment is the numerical 
calculus. We will discuss this in detail in the forthcoming paper. Now 
we are able to employ the truncated moment method in calculation of the 
contribution to moment integrals coming from the region of very small x. 
Using the definition of the truncated moment (|3.1|) we can in a very easy 
way find the double truncated moments of f NS , gi S or -F^ 5 : 

X 2 1 1 

J dxx n - 1 f NS (x,k 2 ) = J dxx n - x f NS {x,k 2 )- J dxx n - l f NS {x,k 2 ) (4.14) 

Xl X\ x 2 

etc. In this way for the partial Bjorken sum rule, according to (|4.5j) . we 
have: 

J dxg^ s (x, k 2 = 10) = JiOn, 1, 10) - h(x 2 , 1, 10) (4.15) 

where x\ and x 2 are both very small: 

xi < x 2 < 10~ 2 (4.16) 
We have thus estimated a contribution: 

AI BSR ( Xl ,x 2 ,k 2 ) = J dxg? s (x,k 2 ) (4.17) 
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and found that AI B sr(W~ a , 1(T 2 , 10) = 0.011, AJ B Sii(10 -5 , HT 3 , 10) = 
0.0013 what is equal to around 5% and 0.6% respectively of the value of the 
full sum (|4.13|) . In our above estimation we have assumed the simple Regge 
input parametrization of g® NS (|4.11|) . 

5. Summary and conclusions. 

In our paper we have derived the integral equation for truncated mo- 
ments of the unintegrated nonsinglet structure function f NS (x,Q 2 ) in the 
case of double logarithmic ln 2 x approximation. Analytical results at fixed 
a s for truncated moments of F^ s and g± have been presented. We have 
received the clear, new solutions and what is important, in the limit xq —* 
our calculations confirm the well known analytical result for full moments 
of f NS (x, Q 2 ). The resummation of ln 2 x terms at small x goes beyond the 
standard LO (and even NLO) QCD evolution of spin dependent parton den- 
sities. In the case of unpolarized structure functions, the Regge behaviour, 
originating from the nonperturbative contribution is the leading one at small 
x. Therefore the double logarithmic approximation is very important par- 
ticularly for the polarized structure function g\ , which at low x is dominated 
just by ln 2 x terms. In our paper we have obtained analytical solutions for 

truncated Mellin moments of g^ s (and F 2 N too) . Dealing with truncated 
l 

moments at xq\ J dxx n ~ l g\{x, Q 2 ) one can avoid uncertainty from the un- 

measurable very small x — > region. In this way the theoretical predictions 
of QCD analysis for structure functions at small x can be compared with 
experimental data without necessity to extrapolate results into nonavailable 
x < xq range. In our analysis for nonsinglet structure functions g^ s and 
F2 we have found their truncated at xq moments, what could be helpful in 
the estimation of Bjorken (g r j V5 ) and Gottfried (F^ 3 ) sum rules. We have 
estimated the contribution to the Bjorken sum rule from the very small x 
region (10 -4 < x < 10~ 2 ) and found it to be around 5% of the value of the 
sum. Evolution equations, we have derived for truncated Mellin moments 
of f NS (x,Q 2 ) generate correctly the leading small x behaviour but do not 
describe Q 2 evolution at large values of x. In the integrals for moments 
of structure functions one need to have values of g% or F2 from the entire 
range of x: xq < x < 1 so for large x one should include in the formalism 
LO Altarelli-Parisi (DGLAP) evolution. The aim of our paper was to focus 
attention on the truncated moments technique within ln 2 x approach itself. 
We found an analytical solution for the truncated moments of structure 
functions within ln 2 x approximation. This is an important stage before 
further investigations. The next step in the improvement of our analysis 
will be including LO DGLAP evolution into equations generating double 
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logarithmic terms ln 2 x for truncated moments of the nonsinglet structure 
functions. This should give a proper determination of truncated and thus 
experimentally verifiable sum rules. 



Acknowledgements 

The work presented here is the result of discussions of one of the authors 
(D.K.) with Jan Kwiecihski. He was, as usual, very helpful, patient and 
stimulating. We are greatly indebted to him for many years of teaching us 
QCD. It is painful to us to express our appreciation of the help we received 
from Jan Kwiecihski after he is gone. 



Appendix A 

Analytical solution of the evolution equation for the full moments of the 
function f NS (x,k 2 ) generating double logarithmic In 2 x effects at small x. 

After double differentiation of (|2.15|) with respect to Ink 2 one obtains: 
*I^= n B J^-.J V) (A,) 

with 

t = lnk 2 (A.2) 
Hence function f NS (n, k 2 ) has a form: 

f NS (n,k 2 )=(^J H(n,a s ) (A.3) 

where H(n,a s ) from initial conditions should be: 

H(n, a,) = f NS (n)R(n, d s ) (A.4) 

and 7 is given by (|2.19[) - (|2.20[) . After inserting the solution l|A.3 |) - (|A.4|) into 
(|2.15|) one finds that R is connected with a s , n and 7 via: 

R(n,ds) = ^- (A.5) 
a s 

Finally, the solutions (|A.3|) - ()A.5|) for n-th moment f NS (n, k 2 ) giving (|2.18j) . 
imply via inverse Mellin transform 

c+ioo 



f»*(x,k z ) = / dnx- n f Nb (n,k 2 ) (A.6) 

2m 



C—lOO 
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the following small x behaviour for f NS (x,k 2 ^- 



f NS (x,r) 



where 



4«-^ (n0lnT + 1)( x r(ln x) 



(A.7) 



(Ai 



and no defined in (|2.2(Jf) denotes branch point singularity of the anomalous 
dimension 7. 

Appendix B 

Analytical solution of the evolution equation for truncated moments of the 
function f NS (x,k 2 ) generating double logarithmic In 2 x effects at small x. 

From the definition of the n-th (n 7^ 0) truncated moment of f(x,k 2 ) 
()3.1|) one can read the following relation: 



df N *(x ,n,k 2 ) _ n df Nb (x Q ,0,k 2 ) 

— %n 



(B.l) 



dx u dx 

On the strength of this above relation we can replace f NS (x , 0, k' 2 ) in (HUH) 



via: 



,v. f , „ ..-„,Y^. . „ / ^/^^ (B.2) 



f NS (x ,0, k' 2 ) = XQ n f (x ,n,k' 2 ) - n 



*0 



and hence we get for ()3.8|) : 



/ 7V5 (x ,n,^) = / JVJ (x ,n) + 



ATS/ 



a, /■ W 2 



// j fc /2 

fc 2 



*1 
k' 2 



f NS (x ,n,k») 



'2\ 



1 



+a s x 



k 2 /xo 



fc' 2 j y 

fc 2 %0 



Now we "guess" the solution of (|3.8|) in the form: 

ATS 



f NS (x ,n,k 2 ) = / JV °(x ,n) [ ^ ) 5(^o,™,a s ) 



fc_ 2 V 

k 2 



(B.3) 



(B.4) 
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where we postulate the same anomalous dimension 7 to be in agreement 
with the result for the full moment: 

f NS (x = 0,n,k 2 )=f NS (n,k 2 ) (B.5) 



After insertion l)B.4Jl into (|B.3|) we can find a simple form of the auxiliary 
function g(xo,n, a s ): 

ft 

g(x ,n,a s ) = ■ (B.6) 

n — 7 + 7Xq 

or equivalently: 

R 

g(%n ' a " s) = i + (fi-iK (R7) 

with R = R(n, a s ) given by (|3.10|> . 
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